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We review the properties of dust in protoplanetary disks around optically visible pre-main 
sequence stars obtained with a variety of observational techniques, from measurements of 
scattered light at visual and infrared wavelengths to mid-infrared spectroscopy and millimeter 
interferometry. A general result is that grains in disks are on average much larger than in the 
diffuse interstellar medium (ISM). In many disks, there is evidence that a large mass of dust 
is in grains with millimeter and centimeter sizes, more similar to "sand and pebbles" than to 
grains. Smaller grains (with micron-sizes) exist closer to the disk surface, which also contains 
much smaller particles, e.g., polycyclic aromatic hydrocarbons. There is some evidence of a 
vertical stratification, with smaller grains closer to the surface. Another difference with ISM is 
the higher fraction of crystalline relative to amorphous silicates found in disk surfaces. There is 
a large scatter in dust properties among different sources, but no evidence of correlation with 
the stellar properties, for samples that include objects from intermediate to solar mass stars 
and brown dwarfs. There is also no apparent correlation with the age of the central object, 
over a range roughly between 1 and 10 Myr. This suggests a scenario where significant grain 
processing may occur very early in the disk evolution, possibly when it is accreting matter from 
the parental molecular core. Further evolution may occur, but not necessarily rapidly, since we 
have evidence that large amounts of grains, from micron to centimeter size, can survive for 
periods as long as 10 Myr. 



1. INTRODUCTION 

Young stars are surrounded by circumstellar disks, made 
of gas and dust. Some of these disks will form planets, and 
one of the open questions of the present time is to under- 
stand which disks will do it, which will not, and why. In 
this chapter, we will refer to all circumstellar disks as pro- 
toplanetary, even if some (many) may only have the poten- 
tial to evolve into a planetary system, and some (many) not 
even that. 

The mass of protoplanetary disks is dominated by gas. 
Still, the solid component (dust grains) is of great impor- 
tance. Observations of light scattered on dust grains and 
their thermal emission remain key diagnostic tools to detect 
disks and to characterize their structure. Dust grains play 
an active role in determining the thermal and geometrical 



structure of disks because their opacity dominates over the 
gas opacity whenever they are present. Furthermore, grains 
shield the disk midplane from energetic radiation, thereby 
influencing the ionization structure, possibly leading to a 
"dead" zone where the magneto-rotational instability can- 
not operate. The formation of dust grains is a phase tran- 
sition which provides solid surfaces important for chemical 
reactions and the freeze-out of molecular components such 
as CO and H2O in the colder parts of disks. Finally, the 
solid mass of the disks is important because the dust grains 
are the building blocks for the formation within the disk of 
planetesimals and eventually planets. 

Grains in protoplanetary disks are very different from 
grains in the diffuse interstellar medium; their properties 
change from object to object and we believe that in each 
object they evolve with time. Without a comprehensive 
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understanding of dust properties and evolution, many disk 
properties cannot be understood. An immediate example 
is the determination of disk solid masses from millimeter 
continuum emission, which requires knowledge of the dust 
absorption coefficient for a specific grain composition and 
structure, and of the grain temperature. 

Protoplanetary disks interest us also because they are 
the place where planets form. Grains are thought to be 
the primary building blocks on the road to planet forma- 
tion, and the continuing and growing interest in this field 
is largely motivated by the need to understand the diversity 
of extrasolar planetary systems and reasons for this diver- 
sity. Earlier reviews by Weidenschilling and Cuzzi (1993) 
and Beckwith et al. (2000) in Protostars and Planets III 
and IV summarize what was known at the time. Since Pro- 
tostars and Planets IV we have seen a lot of progress in 
our knowledge of grain properties in disks, thanks to high- 
resolution data provided by infrared long-baseline interfer- 
ometry, the high sensitivity of Spitzer and ground-based 10 
m-class telescopes, millimetre interferometers such as PdB 
and OVRO, and the VLA capabilities at 7 mm wavelength. 
Data exist now for a large number of disks, around stars 
with very different mass and luminosity, from intermediate- 
mass objects (Herbig AeBe stars, or HAeBe in the follow- 
ing) to T Tauri stars (TTS) and brown dwarfs (BDs). On the 
theoretical side, there is a revived effort in modeling grain 
processing in disks (via coalescence, sedimentation, frag- 
mentation, annealing etc.) and its relation with the chemi- 
cal and dynamical evolution of the gas. Although it is still 
very difficult to integrate observations and theory in a quan- 
titative description of disk evolution, we can expect major 
advances in the near future. 

In this review, we concentrate on observational evidence 
for grain evolution with a special emphasis on grain growth 
and mineralogy. As already mentioned, grain properties 
and their distribution within the disk affect many observable 
quantities. These aspects will be covered in other chapters 
of this book, and we will concentrate on observations that 
measure directly the grain physical and chemical structure. 

Prior to that, we will review very briefly the most im- 
portant processes that control the grain properties (Sec- 
tion 2). We will then outline the various observational tech- 
niques and their Umitations (Section 3), discuss evidence 
for grains of different size, from Polycyclic Aromatic Hy- 
drocarbons (PAHs; Section 4) to micron-size (Section 5) 
and centimeter-size (Section 6). We will discuss grain min- 
eralogy in Section 7. Section 8 will summarize the main 
conclusions one can derive from the observations and out- 
Une some of the open questions we need to address in the 
future. 

2. GRAEV GROWTH: WHY AND HOW 

Grains in the diffuse interstellar medium (ISM) are very 
likely a mixture of silicates and carbons, with a size dis- 
tribution from ~ 100 A to maximum radii of ~ 0.2 — 0.3 
Hm. Smaller particles, most likely PAHs and very small car- 



bonaceous grains, are also present (e.g., Draine, 2003). The 
composition of dust in molecular clouds is similar. Models 
of dust evolution in coUapsing cores (Krugel and Sieben- 
morgen, 1994; Ossenkopf and Henning, 1994; Miyake and 
Nakagawa, 1995) predict only minor changes, as confirmed 
by the observations of Class and Class 1 objects (e.g., 
Beckwith and Sargent, 1991; Bianchi et al, 2003; Kessler- 
Silacci et al., 2005). 

Major changes occur once dust is collected in a cir- 
cumstellar disk, where the pristine interstellar grains may 
grow from sub-micron size or even smaller particles, 
to kilometer- sized bodies (planetesimals), and eventually 
planets. This process is driven by coagulation of smaller 
particles into larger and larger ones, and has a long his- 
tory of theoretical and laboratory studies, which it is not the 
scope of this paper to analyze. We refer for a critical dis- 
cussion to the most recent reviews of Henning et al. (2005) 
and the chapter by Dominik et al. However, the relevance 
of the observations of grain properties we will discuss in 
the following is better understood in the context of grain 
growth and planetesimal formation models, which we will 
therefore briefly summarize here. 

Grain growth to centimeter- and even meter-sizes is 
mainly driven by colhsional aggregation, although gravita- 
tional instabilities may play a role in over-dense dust re- 
gions. The relative velocities leading to such collisions 
are caused by the differential couphng to the gas motion 
(see, e.g., Weidenschilling and Cuzzi, 1993; Beckwith etal, 
2000). This immediately demonstrates that grain growth 
cannot be understood without a better characterization of 
the gas velocity field in disks, especially the degree of tur- 
bulence. For sub-micron-sized grains this is not a real issue 
and their growth by Brownian motion is reasonably well 
understood thanks to numerical simulations and extensive 
laboratory experiments (see, e.g., Blum, 2004 and the chap- 
ter by Dominik et al.). The relative velocities of grains in 
the 1-100 (Um size range are of the order of 10~^ m/s to 
lO"** m/s, low enough for sticking. The outcome of this 
early growth process is a relatively narrow mass distribution 
function with fractal aggregates having open and chain-Uke 
structures. 

While growing in mass, particles also start to sediment 
(due to the vertical component of the stellar gravitational 
field, Schrapler and Henning, 2004, Dullemond and Do- 
minik, 2005 and references therein) and the relative veloc- 
ities may eventually become higher leading to the com- 
paction of aggregates. Above the compaction limit, a run- 
away growth can be expected, where a few large aggre- 
gates grow by colUsions with smaller particles. In this 
regime an exponential increase of mass with time can be ex- 
pected. The coupling between sedimentation, grain growth 
and structure of particles remains the next challenge for pro- 
toplanetary dust models which have to deal with the transi- 
tion from open to compact aggregates. For grains of boulder 
size (larger than 1 m), the relative velocities remain more 
or less constant at about 50 m/s, challenging growth mod- 
els because this velocity is above the destruction threshold 
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(but see Wurm etai, 2005). Furthermore, such boulders are 
rapidly transported to the central star within a few hundred 
years. There are quite a number of proposals to solve the 
problem of the growth barrier, none with complete success 
yet. 

In spite of the many uncertainties, the models converge 
on predicting that grains can grow to very large sizes and 
sediment under the effect of the stellar gravity, leaving a 
population of smaller grains closer to the disk surface and 
of increasingly larger bodies toward the midplane. Grain 
growth should be much slower in regions further away from 
the star than in regions closer to it, which should be com- 
pletely depleted of small grains. When only coagulalion 
mechanisms are considered, small grains are removed very 
quickly, leading to the complete disappearance of the dusty 
disks on timescales much shorter than the age of observed 
disks. Mechanisms that replenish disks of small grains are 
clearly required, for example aggregate fragmentation. Ra- 
dial and meridian circulation can also play an important role 
in grain evolution. 

The process of growth of grains from submicron to km- 
size we have described assumes that the original solid mass 
never returns to the gas phase. This may not be the case, as 
the temperatures in the inner disk are high enough to evap- 
orate grains. If, over time, radial circulations carry close 
to the star most of the disk material, all the original solid 
mass is destroyed and condensed again in the disk. If so, the 
newly formed grains may be very different from those in the 
ISM, even if actual growth does not occur, or is important 
only above a certain (large) size. Gassification and recon- 
densation of grains in the inner disk, coupled with strong 
radial drifts, has been suggested (Gail, 2004) as an expla- 
nation for the large fraction of crystalline grains seen in sev- 
eral disks (see Section 7), and may play an important role 
also in determining the grain size distribution. 

For a long time, the only test to planet formation theories 
was our own Solar system. This has changed completely 
over the last decade, first with the discovery of the variety 
of planetary systems that can actually form, but also with 
a much better knowledge of the grain properties and distri- 
bution in disks which are the progenitors of planetary sys- 
tems. Although our knowledge and understanding of this 
second aspect is still hmited, and the answers to some key 
questions (e.g., which planetary systems will form and from 
which disks) elusive, we will show in the following that the 
observations may already provide important constraints to 
theory. 

3. OBSERVATIONAL TECHNIQUES 

Grain properties in disks can be explored with a variety 
of techniques and over a large range of wavelengths, from 
visual to centimeter. Before reviewing the results and their 
implication, it is useful to sunnmarize briefly the capabihties 
and limitations of the different techniques. 

The most severe limits to the characterization of grains 
in disks do not come from observational limits, such as spa- 



tial resolution or sensitivity (although these too can be a 
problem), but from the physical structure of disks. Their 
extremely high optical depth allows us to get information 
only on grains located in a r = 1 layer, which at optical and 
IR wavelengths contains a tiny fraction of the dust mass. 
Moreover, if one measures dust emission, the observed flux 
is strongly biased by the temperature dependence of the 
Planck function. 

The only way to measure the properties of the bulk of 
the dust mass is to go to longer wavelengths, where an in- 
creasing fraction of the disk becomes optically thin. Inter- 
ferometric observations at millimeter and centimeter wave- 
lengths have provided the strongest evidence so far that 
most of the original solid mass in protoplanetary disks has 
grown to cm- size by the time the central star becomes opti- 
cally visible. However, one should remember that this tech- 
nique is limited today by the sensitivity and resolution of 
existing milhmeter interferometers. For example, studies 
of the radial dependence of the dust properties are still very 
difficult. Furthermore, even at these very long wavelengths, 
the regions of the disk closer to the star (typically, up to few 
AUs) are optically thick and thus this very interesting iimer 
portion of the midplane remains inaccessible at aU wave- 
lengths. 

The mid-infrared spectral region, roughly between 3 and 
100 /im, is rich in vibrational resonances of abundant dust 
species, sihcates in particular. The spectral region between 
3 and 20 /zm contains also prominent R\H C-C and C-H 
resonances. The wavelength, spectral shape and strength of 
the resonances are sensitive probes of the chemical com- 
position, lattice structure, size and shape of the particles. 
In disks around optically visible stars these dust features 
are generally seen in emission, having origin in the opti- 
cally thin surface layers, heated by the stellar radiation to 
temperatures higher than the disk midplane (e.g., Calvet et 
ah, 1991; Chiang and Goldreich, 1991 \ Menshchikov and 
Henning, 1997). Infrared spectroscopy roughly probes dust 
particles with sizes of up to a few microns (depending on 
material and wavelength), in the temperature range between 
1500 and about 50 K. In protoplanetary disks, this implies 
that such spectra are sensitive to dust in the innermost re- 
gions, and in the surface layers of the disk only. No infor- 
mation is obtained about chemical and allotropic properties 
of large and/or cold dust grains. 

The 10 jim spectral region deserves special attention be- 
cause it is accessible from the ground, is available for high 
spatial resolution interferometric studies, and a fairly large 
body of data has been collected, covering a wide range of 
stellar mass and age. It contains the strongest resonances of 
both amorphous and crystalline silicates. The silicate emis- 
sion bands have important but somewhat hmited diagnostic 
value. First, the 10 jjm spectral region probes grains of a 
certain temperature range, mostly between 200 and 600 K 
Second, both the amorphous and crystalline sihcate reso- 
nances near 10 /itm only probe presence of grains smaller 
than a few microns in size, the exact value depending on 
the dust species. This limits the diagnostic value of the 10 
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nm region to warm, small silicates mostly located on the 
disk surface in regions within ^10 AU for protoplanetary 
disks around HAe stars, ~ 0.5 - 1 AU for TTS, and ~ 0.1 
AU for BDs. 

At shorter wavelengths (visual and near-IR), scattering 
becomes important and one can observe scattered hght ex- 
tending to very large distances from the star (e.g., Mc- 
Caughrean et ai, 2000). Multi-wavelength images in scat- 
tered light can be used to derive grain properties over a large 
range of disk radii, but hmited to a very narrow layer on the 
disk surface. 

Interesting information on grain properties deeper in the 
disk can be obtained from edge-on disks (silhouette), seen 
as a dark lane against a luminous background. Informa- 
tion on dust properties can be obtained in favorable cases 
by studying the wavelength dependence of the dark lane 
absorption against the background of the scattered light 
at visual and near-infrared wavelengths. Silhouette disks 
should have deep silicate features at 10 /xm in absorption, 
which, contrary to the features in emission, sample the 
colder grains in the outer disk. Sihcate absorption is ob- 
served in few cases, but it is often difficult to exclude con- 
tamination from grains in the surrounding material. 

Multiwavelength studies are obviously necessary to 
completely characterize the grain population, in size, com- 
position and physical structure. So far, this cannot be done 
in disks, where observations at different wavelengths sam- 
ple different regions of the disk. As we will describe in 
the following, there is clear evidence that grains of all sizes 
(from PAHs to cm-size bodies at least) are present in disks. 
However, the information is always partial, and it is very 
difficult to get a global picture, for example, to measure the 
total mass fraction of grains of different size. 

An additional caveat is that none of the techniques we 
have mentioned is sensitive to bodies larger than few cen- 
timeters. Kilometer- size planetesimals can be detected 
through the dynamical perturbations they create, but the de- 
tection of meter-size bodies is practically impossible. 

4. THE SMALLEST PARTICLES: PAHS 

Emission from transiently heated very small particles (at 
3.3, 6.2, 7.7-7.9, 8.6, 11.3 and 12.7 fim etc.) has been 
detected in many disk systems, mostly HAeBe stars (e.g., 
Acke and van den Ancker, 2004) and generally attributed to 
PAHs. For a long time, it has not been clear if the observed 
emission is from PAHs in the disks or from the reflection 
nebulosities associated with many HAeBe. There is now 
convincing evidence in favour of the disk origin. 

Peeters et al. (2002), van Diedenhoven et al. (2004) and 
Sloan etal. (2005) have shown that in HAeBe stars the PAH 
features differ in shape and wavelength when compared to 
the ISM. Meeus et al. (2001) and Acke and van den An- 
cker (2004) found that the strength of the PAH bands corre- 
lates with the shape of the spectrum in the 10-60 /im region. 
Flared disks (i.e., with flux increasing with wavelength 
in the 10-60 fiia range; group I in the Meeus et al. (2001) 



classification) have strong PAHs features, while flat disks 
(i.e., with Fi, decreasing with wavelength in the 10-60 /xm 
range; group II) have no or very weak emission. This trend 
has been analyzed by Habart et al. (2004a), who have com- 
puted disk models which include PAHs and shown that the 
observed correlation between the strength of the PAH bands 
and the shape of the SED is well explained in terms of the 
solid angle the disk subtends as seen from the star. In a 
few stars, other features, at 3.43, 3.53 /xm, from transiently 
heated carbonaceous materials, possibly identified as nano- 
diamonds (Guillois et al, 1999), have also been shown to 
have a disk origin {Van Kerckhoven et al., 2002). 

Recently, PAH and nanodiamond emission has been spa- 
tially resolved in several objects (van Boekel et al, 2004b; 
Ressler and Barsony, 2003; Habart et al, 2004b, 2006); 
the emitting region has a size consistent with disk model 
predictions (Habart et al, 2004b, 2006; see Fig. 1). The 
emission of such transiently heated species (which are "hot" 
whenever excited) is much more extended than that of the 
adjacent continuum, emitted by larger grains in thermal 
equihbrium with the radiation field which become rapidly 
cool as the distance from the star increases. Using images 
in the features of PAHs and nanodiamonds it has been pos- 
sible to measure disk inclination and position angle on the 
sky {Habart et al, 2004b) at spatial scales of few tens AU. 

The presence of very small particles on disk surfaces has 
a strong impact on the gas physical properties, since they 
contribute a large fraction of the gas heating via the photo- 
electric effect and may dominate the H2 formation on grain 
surfaces. PAHs are thus an essential ingredient of disk gas 
models (e.g. chapter by Dullemond et al). However, the 
observations have severe limitations, as they can only probe 
matter at large altitude above the disk midplane, so that we 
do not know if PAHs survive deeper into the disk; moreover, 
the intensity in the features decreases with radius , and is 
well below current observational capabihties for distances 
>40 AU (at least for the 3.3 jim feature). 

There are also limitations in the current models, which 
need to explore a much wider range of parameters, both in 
disk and PAH properties. 

5. FROM SUB-MICRON TO MICRON GRAINS 

5.1 Scattering and Polarization 

Observations at visible and near-infrared wavelengths of 
the dark disk silhouette seen against the background light 
scattered by the disk surface provide direct view of disks 
around young stars. The properties of the scattered radia- 
tion and the fraction of polarization can be used in principle 
to measure the size of the scattering grains. However, we 
know from laboratory work that particle sizing based on the 
angular distribution of scattered light or the measurement 
of polarization is not without problems and usually only 
works well for monodisperse distributions of spheres. In as- 
tronomical sources, the analysis of scattered radiation and 
polarization is even more complicated because it is not pos- 
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Fig. 1. — Intensity profile of the 3.3 /im PAH feature in 
the HAe star HD 169142, observed with NAOS/CONICA 
on the VLT (Habart et al, 2006). The top panel shows the 
intensity in the feature after continuum subtraction (thick 
solid line) and in the adjacent continuum (thin solid line); 
the dashed line plots the point-spread function as measured 
on a nearby unresolved star. The PAH feature is clearly 
spatially resolved, with FWHM size of 0.3 arcsec (roughly 
40 AU at the distance of the star). The continuum is unre- 
solved. The bottom panel compares the observations with 
disk model predictions for the feature (dot-dashed line) and 
the continuum (dotted line). Details of observations and 
calculations, including model parameters, are in Habart et 
al. (2004a, 2006). 

sible to obtain measurements for all scattering angles. In 
addition, polarization measurements are often technically 
challenging which may explain why we only have limited 
data using this technique and grain size information from 
such observations is still scarce. 

Scattering is well described by the 4x4 scattering ma- 
trix (Mueller matrix) transforming the original set of Sto- 
kes parameters into a new set after the scattering event. The 
matrix elements are angular-dependent functions of wave- 
length, particle size, shape, and material. For particles 
that are small compared with the wavelengths, the scattered 
light is partially polarized with a typical bell-shaped angu- 
lar dependence with 100% polarization at a scattering angle 
of 90 degrees. The polarization degree depends only on the 
scattering angle and not on particle size. This practically 
means that the radius of small spherical particles cannot be 
determined from scattering measurements. 

For very small particles scattering does not vary a lot 
with direction. This changes drastically for larger spheres 
with more ripples and peaks in the scattering pattern, for 



which forward scattering becomes very important. In ad- 
dition, the peak polarization decreases and moves to larger 
scattering angles. The very structured scattering pattern for 
large spheres is used for experimental particle sizing, but 
the characteristic ripples disappear if a grain size distribu- 
tion is considered instead of a single size population. For 
very large spheres we again find smooth curves. Further- 
more, non-spherical particles may behave very differently. 
We should also note that an ensemble of non-identical scat- 
tering particles leads to depolarization of scattered light, 
adding to the complexity of our understanding of the polar- 
ization degree (in terms of grains sizes) of radiation coming 
from a disk. We refer to Voshchinnikov and Kriigel (1999) 
for Mie calculations, demonstrating the effect of different 
grain sizes. 

One may conclude that solving the inverse problem of 
determining grain sizes from scattering and polarization ob- 
servations is hopeless. This is certainly not true. A grow- 
ing number of objects are studied using multi-wavelength 
imaging and polarization techniques. For our purposes, it 
is sufficient to note that, in all cases, the evidence points 
toward grain growth on the disk surface to sizes of up to 
a few microns; often, there is also evidence of sedimenta- 
tion, i.e., that larger grains are closer to the disk midplane. 
Among the more recent studies, we note the work of Lucas 
et al. (2004), who nicely demonstrated the power of a de- 
tailed analysis of high-resolution imaging polarization data 
for the case of HL Tau. They found that silicate core-ice 
mantle grains with the largest particles having radii slightly 
in excess of 1 /im best fit the data. A similar study has 
been performed for the ring around GG Tau by Duchene et 
al. (2004). They found again a slight increase in the grain 
size towards micron-sized grains, and evidence for a verti- 
cal stratification of dust. The surface layers, located 50 AU 
above the ring midplane, contain dust grains that are con- 
sistent with being as small as in the ISM, while the region 
of the ring located 25 AU from the midplane contains sig- 
nificantly larger grains (> 1 /im). This stratified structure is 
likely the result of vertical dust settling and/or preferential 
grain growth in the densest parts of the ring. The first scat- 
tering measurements at mid-infrared wavelengths obtained 
for the T Tauri star HK Tau B (McCabe et al, 2003) showed 
similar results. 

Evidence of growth of grains to sizes of up to a few mi- 
crons and of a vertical stratification of grain sizes within 
the disk surface has also been inferred in silhouette disks 
in Taurus {D'Alessio et al., 2001) and Orion {Throop et 
al., 2001; Shaping et al., 2003; see also McCaughrean et 
al., 2000 and references therein), by studying the behaviour 
with wavelength of the translucent edges of the dark lane. 

5.2 Mid-Infrared Spectroscopy 

High quality mid-IR spectra are available for an increas- 
ing number of objects. First ISO and now Spitzer are pro- 
viding us with data of a quality that could not be achieved 
from the ground, even if the spatial resolution is in both 
cases too low to resolve the emission. Spitzer spectra are 
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becoming available in increasing number at the time we 
write, and we expect further improvement in the statistics 
and quality of the data, for low mass objects in particular. 
In the following, we provide a summary of the silicate prop- 
erties, as derived from the available observations. 

Most of the information is obtained from observations of 
the 10 fim emission feature. The spectra show a large vari- 
ation in the strength and shape of this feature, for objects of 
all mass, from HAe stars to TTS and BDs. Fig. 2 plots a 
representative selection of HAe stars, TTS, and BDs, which 
shows the variety of observed profiles. In some cases, the 
shape of the feature is strongly peaked at about 9.8 fxm, as 
for small 1 /im) amorphous grains in the ISM. In other 
objects, the feature is very weak with respect to the contin- 
uum, and much less peaked. In other cases, many narrow 
features, typical of crystalline siUcates, are clearly visible, 
overimposed on the smoother and broader amorphous sil- 
icate emission. These features will be discussed in detail 
in Section 7. In spite of these large variations, there is a 
general trend, valid for objects of all mass, first identified 
by van Boekel et al. (2005), who noted that the shape and 
strength of the silicate band are correlated, so that weaker 
features are also flatter This is consistent with the growth 
of siUcate grains from sub-micron to micron sizes. If grains 
grow further, the siUcate emission disappears; there are in 
fact a few objects which show no silicate emission, most 
likely because all small sihcates have been removed (e.g. 
Meeus et al, 2003). 

A useful way to summarize the properties of the 10 /im 
feature of a sample of objects is shown in Fig. 3. For each 
object, the shape of the 10 /itm silicate emission feature is 
characterized by two parameters, the ratio of the flux at 
the peak of the feature over the continuum (FpeoA: / Fcont) 
and the ratio of the continuum subtracted 11.3 /xm over 
the 9.8 /im flux. Amorphous silicates of increasing size 
have weaker and flatter features (cf. Fig. 2), with values 
of Fpeak /Fcont decreasing from 3 — 4 (for sizes <Cl /xm) to 
^ 1, for grains larger than few microns. At the same time, 
the ratio -Fn a/i^g g increases from ~ 0.5 — 0.6 to ~ 1. 
Larger values of the ratio -F11.3/-F9.8 are due to a significant 
contribution from forsterite (i.e., crystalline Mg2Si04). 

Note that the ratio Fp^ak/ Fcont depends not only on 
grain cross section but also on other quantities (e.g., the 
disk inclination and geometry) which affect differently fea- 
ture and continuum emission (e.g., Chiang and Goldreich, 
1999). Also, changes in grain properties other than size 
(e.g., porosity; see Voshchinnikov et al, 2006, Min et al, 
2006) can reproduce part of the observed trend, and make 
the simple analysis in terms of grain size more uncertain. 
However, the interpretation of the 10 /im feature profiles in 
terms of growth of the grains in the disk surface is convinc- 
ing, and in agreement with the results from other kinds of 
observations discussed in Section 5.1. 

Additional information will be derived in the future from 
the properties of the weaker 20 /xm silicate band, which 
samples slightly larger grains, and regions on the disk sur- 
face further away from the star. Spitzer spectra are becom- 



ing available for a rapidly growing number of objects, and 
the first results seem to confirm the results obtained from 
the 10 /xm feature analysis {Kessler-Silacci et al, 2006; 
Bouwman et al, 2006). 

6. GROWTH TO CM-SIZE GRAINS 

As discussed in Section 3, at (sub-)millimeter wave- 
lengths and beyond the dust emission in disks begins to be 
moderately optically thin. At these wavelengths it is thus 
possible to probe the bulk of the dust particles, which are 
concentrated in the disk midplane. The spectral energy dis- 
tribution of the continuum emission can be directly related 
to the dust emissivity, which in turn is related to the grain 
properties. The variation of the dust opacity coefficient per 
unit mass with frequency, in the millimeter wave range, can 
be approximated to relatively good accuracy by a power 
law: k ~ A"''. The value of the /? exponent is directly 
related to the dust properties, in particular to the grain size 
distribution and upper size cutoff; in the hmit of dust com- 
posed only of grains much larger than the observing wave- 
length, the opacity becomes grey and (3 = 0. 

It has been known for many years that pre-main se- 
quence disks have (sub-)millimeter spectral indices a 
(Fu oc A^") shallower than prestellar cores and young pro- 
tostars (Beckwith and Sargent, 1991). If the disk emission 
is optically thin and hu/kT^ <ti 1, the observed values of 
a would translate very simply in opacity power law indices 
/3 = a — 2; with observed typical values of a < 3, this 
would imply /3 ~ 1 or lower. This immediately suggests 
that disk grains are very different from ISM grains, which 
have P 1.7 (Weingartner and Draine, 2001). 

This result, however, has been viewed with great caution, 
since the effect of large optical depth at the frequency of the 
observations cannot be ruled out for spatially unresolved 
observations. The extreme case of an optically thick disk 
made of ISM grains has a = 2, the same value as an opti- 
cally thin disk with very large grains, for which (3 = (see, 
e.g., Beckwith and Sargent, 1991; Beckwith et al, 2000). 
The combination of optical depth and grain properties can 
be disentangled if the emission is resolved spatially and the 
disk size at the observed wavelength is measured (see, e.g., 
TesfiefaL, 2001, 2003). 

There are now several disks whose millimeter emission 
has been imaged with interferometers. A major step for- 
ward has come from the use of 7 mm VLA data, which, 
in combination with results at shorter wavelengths (1.3 and 
2.6 mm) with PdB and OVRO, provides not only high spa- 
tial resolution and larger wavelength range (minimizing the 
uncertainties on a), but also the possibility of probing pop- 
ulation of grains as large as a few centimeters. Addition- 
ally, the VLA at centimeter wavelengths has been useful to 
check whether the dust emission at shorter wavelength is 
contaminated by free-free emission (e.g. Testi et al, 2001; 
Wilner et al. 2000; Rodmann et al. 2006). We summarize in 
Fig. 4 the results for the objects with 7 mm measurements 
from Natta et al (2004) and references therein and Rod- 
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Fig. 2. — Observed profiles of the 10 /im silicate feature. The first three Panel (from left to right) show a selection of 
profiles for HAe stars (van Boekel et ai, 2005), TTS (Bouwman et ai, 2006), and BDs in Chamaleon (Apai et ai, 2005). 
All the profiles have been normalized to the 8 /im flux and shifted for an easier display. In each Panel, they have been 
ordered according to the Fpeak/Fcont value. The dashed vertical lines show the location (at 9.8 /im) of the peak of small 
amorphous olivine (Mg2Si04) and of the strongest crystalline feature of forsterite at 11.3 /im. Amorphous silicate grains 
of composition Mg2a;Fe2_2a;Si04, generally with a; <C 1 have the same composition of olivines, and are are often referred 
to as "amorphous olivines". Although this is not correct (olivine is crystalline), we will sometime use this definition. 
The Panel to the extreme right shows for comparison a selection of profiles of laboratory silicates. Starting from the top, 
amorphous olivine with radius 0.1, 1.5, 3 and 5 /im, as labelled. In all cases, we have added the same continuum (30% 
of the peak of the 0. 1 /im amorphous olivine), and normalized as the observed profiles. The two bottom curves show the 
profile of small crystalline grains, forsterite (Mg2Si04) and enstatite (MgSiOs), respectively; larger crystalline grains have 
broader and weaker features (e.g., van Boekel et ai, 2005). Several of the features displayed by crystalline sihcates can be 
identified in some of the observed spectra. 



mann et al. (2006). For these objects, the contribution from 
gas emission has been subtracted from the millimeter emis- 
sion, and the spectral index a has been computed over the 
wavelength interval 1.3-7 mm; in three cases the 7 mm flux 
is an upper limit, once the gas emission has been subtracted 
from the total. The values of (3 have been derived fitting 
disk models to the data. Note that these model-fitted j3 val- 
ues are slightly larger than the optically thin determination 
(3 — a — 2, typically by about 0.2. This difference is due to 
the fact that the emission of the inner optically thick disk, 
albeit small, is not entirely negligible. The uncertainty on 
the model-derived (3 is of about ±0.2 {Natta et al, 2004; 
Rodmann et ai, 2006). The figure shows also the results 
for those objects for which no resolved maps exist, but have 
measurements of the integrated flux at 7 mm. For these, 
/3 as been derived from the measured a assuming that the 
disk is larger than ^100 AU, sufficient to make the optically 
thick contribution small. 

If we consider only the 14 resolved disks, which have all 
radii larger than 100 AU, (3 ranges from 1.6 (i.e., very sim- 
ilar to the ISM value), to 0.5. Ten objects have /3 < 1, i.e.. 



much flatter than ISM grains. The unresolved disks behave 
in a similar way; the extreme case is UX Ori, with /3 ^ 0. 
Unfortunately, the UX Ori disk has not been resolved so far; 
at the distance of ^ 450 pc, and, with an integrated flux at 
7 mm of only 0.8 mJy (Testi et ai, 2001), it remains a tan- 
talizing object. High and low values of f3 are found both for 
HAe stars and for TTS, and there is no apparent correlation 
of [3 with the stellar luminosity or mass. Similarly, we do 
not observe any correlation of (3 with the age of the star. 
This last quantity, however, is often very uncertain, and the 
observed sample limited. 

Once we have established that the grains originally in the 
collapsing cores, for which /3 ~ 1.7 — 2, have gone through 
a large degree of processing, one wants to derive from /3 the 
properties of the actual grain population, and in particular 
its size. The properties of the millimeter opacity of grains 
of increasing size have been discussed by, e.g., Beckwith 
et al. (2000) and references therein. Fig. 5 illustrates the 
behaviour of (3 computed between 1 and 7 mm for a popu- 
lation of grains with size distribution n(a) cx a^'^ between 
amin and amax- The index (3 is plotted as function of amax 
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Fig. 3. — The ratio of the continuum-subtracted flux at 11.3 
fim over that at 9.8 /im (i^n.a/Fg.g) is plotted as function 
of flux at the peak of the 10 /im feature over the contin- 
uum, Fpeak/Fcont- Open and filled symbols show the ob- 
servations for HAe stars (van Boekel et ai, 2005; filled 
dots), TTS (Kessler-Silacci et al., 2006, Bouwman et al, 
2006, Przygodda et al., 2003; squares), BDs (Apai et al, 
2005; stars). All the data but those of Kessler-Silacci et al. 
(2006) and Bouwman et al. (2006) have been re-analyzed 
by Apai et al. (2005) in a homogeneous way. Error bars 
are shown when available in the literature. The solid curve 
shows the result for a dust model which includes amorphous 
olivine of two sizes (1 and 10 /im) in proportion varying 
from all 1 fim (extreme right) to all 10 pm (extreme left), 
and an additional contribution of a 2% mass fraction of 0.1 
/im forsterite grains (from J. Bouwman, private communi- 
cation); the same continuum (25% of the peak flux of the 
smallest grains) is added to all model specti^a. 
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Fig. 4. — The figure plots the values of the millimeter opac- 
ity index (3 (k (x X~^) as function of the stellar luminos- 
ity (top panel) and of the stellar age (bottom panel); data 
fmm Natta et al. (2004) and Rodmann et al. (2006). Filled 
squares are objects where the disk millimeter emission has 
been spatially resolved, open squares objects which have 
not been observed with high spatial resolution, or that were 
found to be not resolved with the VLA (resolution ^ 0.5 
arc sec). 



for different values of q; in all cases, amin ^ 1 mm. One 
can see that, for all values of q, as a^nax increases f3 is first 
constant at the value typical of grains of size -C 1 mm, it 
has a strong and rather broad peak at a„iax ~ 1 mm and 
decreases below the initial value for a„iax^ few mm. How- 
ever, only for q < 3 /3 goes to zero for large a„iax', for 
q > 3, the small grains always contribute to the opacity, so 
that (3 reaches an asymptotic value that depend on q and on 
the f3 of the small grains. For q = 4, the asymptotic value 
is practically that of the small grains. 

The results shown in Fig. 5 have been computed for 
compact segregated spheres of olivine, organic materials 
and water ice (Pollack et al., 1994). The exact values of 
[3 depend on grain properties, such as their chemical com- 
position, geometrical structure and temperature (e.g., Hen- 
ning et al, 1995); examples for different grain models can 



be found, e.g., in Miyake and Nakagawa (1993), Krilgel 
and Siebenmorgen (1994), Calvet et al. (2002), Natta et 
al. (2004). However, these differences do not undermine 
the general conclusion that the observed spectral indices re- 
quire that grains have grown to sizes much larger than the 
observing wavelengths (e.g., Beckwith et al, 2000 and ref- 
erences therein, Draine, 2006). To account for the obser- 
vations which extend to 7 mm, objects with /3 ^ 0.5 — 0.6 
need a distribution of grain sizes with maximum radii of 
few centimeter, at least. The minimum grain radii are not 
constrained, and could be as small as in the ISM; however, 
the largest particles need to contribute significantly to the 
average opacity, which implies that the grain size distribu- 
tion cannot be steeper than n{a) oc Qr-(3^3-5); in any case, 
the large grains contain most of the solid mass. A detailed 
analysis can be found in e.g., Natta and Testi (2004). 
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Fig. 5. — Opacity index (3 for grains with a size distribu- 
tion n{a) oc a^^ between and a„iax- The index (3 
is computed between 1 and 7 mm for compact segregated 
spheres of ohvine, organic materials, water ice {Pollack et 
ai, 1994), and plotted as function of Umax- Different curves 
correspond to different values of q, as labelled. In all cases, 
amin ^ 1 mm. The dotted curve shows /3 for grains with 
single size amax- The histogram on the right shows the dis- 
tribution of the values of /3 derived from millimeter obser- 
vations (see Fig. 4). About 60% of the objects have /3 < 1, 
and a„iax^ 1 cm. Different dust models give results which 
are not very different (e.g., Natta and Testi, 2004). 



Of the objects observed so far, TW Hya is a particularly 
interesting case. Millimeter observations can be fitted well 
with models that include grain growth up to 1 cm (Calvet et 
ai, 2002; Natta and Testi, 2004; Qi et ai, 2004). However, 
this is probably an underestimate of the maximum grain 
sizes in this disk. Wilner et al. (2005) have demonstrated 
that the 3.5 cm emission of this object, contrary to expec- 
tations, is not dominated by gas emission, but by thermal 
emission from dust grains in the disk, and that a popula- 
tion of particles as large as several centimeters residing in 
the disk is needed to explain both the long wavelength SED 
and the spatial distribution of the 3.5 cm emission. 

The result that dust particles have grown to considerable 
sizes, in fact to "pebbles", is thus a solid one and verified 
in many systems. So far, it has been derived as a "global" 
property of the dust grain population of the disk. In order 
to constrain dust growth models, the next logical step is to 
check for variation of dust properties as a function of radius. 
A first attempt at this type of study is shown in Fig. 6. High 
angular resolution and high signal to noise 7 mm (VLA) 
and 1.3 mm (PdBI) maps of the HD 163296 system have 
been combined to reconstruct the millimeter spectral index 



profile as a function of radius. The result suggests that the 
spectral index varies as a function of radius from a ~ 2.5 
in the inner ~60 AU to a ~ 3.0 in the far outer disk. If 
interpreted as differential grain growth, this result implies 
that the outer disk contains less evolved particles than the 
inner disk. 
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Fig. 6. — The HD163296 disk. In the top panel we show 
the IRAM-PdBI contour map of the 1 .3 mm emission from 
the disk, and the approximate shape and size of the scatte- 
red light disk as detected with HST/STIS. The images have 
been rotated so that the major axis of the disk is aligned 
with the abscissa of the plot. The spatial resolution of the 
millimeter observations is shown by the ellipse in the lower 
right box of the upper panel. In the bottom panel the dia- 
monds with error bars show the variation along the disk ma- 
jor axis of the continuum spectral index measured between 
1.3 and 7 mm. In the same panel, the solid line shows the 
values of the spectral index when a correction for the maxi- 
mum possible gas contribution is applied to the 7 mm obser- 
vations. The data suggest a variation of the spectral index as 
a function of disk radius. If interpreted as a variation of the 
grain properties, they indicate that larger grains are found in 
the inner (r < 60 AU) regions of the disk, while the grains 
in the outer disk may still be in a less evolved stage. 



7. MINERALOGY 

The Infrared Space Observatory (ISO) and the latest 
generation of mid-infrared spectrographs on large ground- 
based telescopes have revealed the presence of many 
new dust species, in addition to amorphous silicates, in 
the protoplanetary disks surrounding pre-main-sequence 
stars. The position and strength of the observed emission 
bands match well those of Mg-rich and Fe-poor crystalline 
silicates of the olivine (Mg22;Fe2-2a;Si04 and pyroxene 
Mga;Fei_a:Si03 families, in particular their Mg-rich end 
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members forsterite and enstatite respectively (x=l). Apart 
from these components, evidence for FeS and Si02 was also 
found {Keller et al, 2002). Spitzer is now providing similar 
results for stars of lower mass, TTS and BDs. 

The presence of these dust species in protoplanetary 
disks impUes substantial processing (chemical and physi- 
cal) of the dust because their abundance, relative to that of 
amorphous silicates, is far above limits set for these species 
in the interstellar medium (e.g., Kemper et al., 2004). 
Therefore, these dust species must be formed sometime dur- 
ing the collapse of the molecular cloud core or in the ac- 
cretion disk surrounding the young star. There are several 
possible mechanisms that may be responsible for grain pro- 
cessing, and it is quite likely that several contribute. In the 
innermost disk regions, where the dust temperature is above 
1000 K, heat will induce a change of the grain lattice order- 
ing (long-range ordering) from chaotic to regular (thermal 
annealing), leading to the transformation of amorphous sil- 
icates into crystalline ones (olivines and pyroxenes). Above 
about 1200-1300 K, chemical equilibrium processes will 
lead to vapourisation and gas-phase condensation of sili- 
cates, mostly in the form of crystalline forsterite. Local 
high-energy processes (shocks, lightning) as well as radial 
mixing may increase substantially the amount of crystalline 
silicates in the cool outer regions of the disk. Note that ther- 
mal annealing of amorphous, presumably Fe-rich silicates 
would lead to Fe-rich crystalhne silicates. Therefore the 
Fe content of crystalhne silicates may be used as an indica- 
tion of the chemical processing that has occurred (gas-phase 
condensation, anneahng, local processes). In this context, 
we stress the importance of improving the still poorly con- 
strained stochiometry of interstellar amorphous silicates, in 
particular the Fe/Mg and Mg/Si ratio. 

The mineral composition of dust in protoplanetary disks 
can be compared directly to that of solar system comets, as- 
teroids and interplanetary dust particles (IDPs). Only one 
high-quahty 2-200 ixm spectrum of a solar system comet, 
Hale-Bopp, is available from ISO. The crystalline silicates 
in Hale-Bopp are very Mg-rich and Fe-poor, similar to 
those in protoplanetary disks. Laboratory studies of IDPs 
of cometary origin also reveal the presence of Mg-rich crys- 
talline silicates. The situation for asteroids is much more di- 
verse. Apart from Mg-rich crystalhne silicates, many aster- 
oids show Fe-containing silicates, probably related to par- 
ent body processing or to nebular processing in the inner 
solar nebula. 

7.1 Spatially unresolved 10 jjxa spectra 

Information on the abundance of crystalline sihcates in 

a large sample of objects of different mass are derived from 
the shape of the 10 /xm feature (mostly through the 1 1.3 /zm 
forsterite peak) in spatially unresolved observations. Fig. 2 
shows how this feature can be prominent in stars of all mass, 
from HAeBe to brown dwarfs. Fig. 3 shows how in many 
objects the flux at 11 .3 fim is stronger than at 9.8 fjm, where 
the emission of small amorphous silicates roughly peaks. A 
strong 1 1.3 /um peak is typical of profiles with a large com- 



ponent of crystalline sihcates, in addition to the amorphous 
ones. This interpretation has been confirmed by the pres- 
ence of crystalhne features at longer wavelengths, when the 
data are available. Contamination from the 11.3 /zm PAH 
feature has been ruled out in many cases, from the absence 
of the other PAH features at shorter wavelengths (see, for 
example, Acke and van derAncker, 2004). 

A determination of the fractional mass abundance of the 
crystalline silicates has been derived by various authors fit- 
ling the observed profiles with mixtures of grains of dif- 
ferent composition (typically, olivines and pyroxenes), size 
(from submicron to few microns) and aUotropic state (amor- 
phous and crystalline). In general, the models assume a 
single temperature for all grains, and that the emission is 
optically thin. These assumptions, and the limited number 
of dust components that are included in models, make the 
estimates very uncertain. However, the results are interest- 
ing. For HAeBe stars, van Boekel et al. (2005) find that 
crystalline silicates are abundant only when the amorphous 
silicate grain population is dominated by large (> 1.5 /im) 
grains. In Fig. 3, there are no points with 11.3/9.8 ratio 
>1 and large values of Fp^ak / Fcont- The fraction of grains 
which is crystalline ranges from 5% (detection limit) to k, 
30%. Note that grains larger than few microns are not ac- 
counted for; if crystalline sihcates are mostly sub-micron in 
size, while the amorphous species have a much broader dis- 
tribution, the degree of crystalhnity can be much lower than 
current estimates. This seems to be the case in Hale-Bopp 
(Min et al, 2005), and may happen in protostellar disks as 
well. In addition, the derived crystalline/amorphous abun- 
dances are spatially averaged values; there are hkely strong 
radial gradients in crystalhnity (see below). As more data 
for lower mass objects become available, we see that they 
show a similar trend, although relatively high crystallinity 
and relatively few large amorphous grains have been seen 
in one very low mass TTS (e.g., Sargent et al, 2005). Re- 
cent Spitzer data for brown dwarfs in Chal show that also 
objects of very low mass have large fraction of crystalhne 
sihcates, from 9 to about 50% {Apai et al, 2005). 

As far as the chemical composition of the silicates is con- 
cerned, there are some general trends, showing that crys- 
talhne pyroxene grains can be found in disks with relatively 
large fractions of forsterite and sihca, but they are typically 
not found in disks with small fractions of forsterite and sil- 
ica {van Boekel et al., 2005; Sargent et al, 2006). 

7.2 Spatially resolved spectroscopy 

The advent of spectrally resolved mid-infrared interfer- 
ometry using large baselines at the Very Large Telescope 
Interferometer (VLTI) has made it possible for the first time 
to study the nature of the dust grains on AU spatial scales. 
The MIDI instrument at the VLTI has been used to study 
the innermost regions of HAeBe star disks, van Boekel 
et al. (2004a; see Fig. 7) show that the inner 1-2 AU of 
the disk surface surrounding three HAe stars is highly crys- 
talline, with between 50 and 100 per cent of the small sil- 
icates being crystalline. However, crystalhne sihcates are 
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also present at larger distance, with a wide range of abun- 
dances. A similar behaviour is found in other objects ob- 
served with MIDI (Leinert et ai, 2006). Clearly, large star- 
to-star variations in the amount and distribution of the crys- 
talline siUcates exist, but it seems that crystalline grains are 
relatively more prominent in the inner than in the outer disk. 
Note that, because silicate emission from cold grains be- 
come very weak, in this case "outer" refers to warm regions 
of the disk between 2 and 20 AU at most (van Boekel et al. , 
2004a). 
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Fig. 7. — Comparison of the siUcate profiles from the iimer 
(1-2 AU) and outer (2-20 AU) disk in three HAe stars (van 
Boekel et al., 2004a). The data have been obtained with the 
mid-IR interferometric instrument MIDI on the VLTI. The 
comparison shows that, for all the three stars, the crystalline 
features are much stronger in the iimer than in the outer 
disk. 

In the star with the highest crystallinity, HD 142527, 
the observations suggest a decrease in the ratio between 
forsterite and enstatite with increasing distance from the 
star. Such a trend in the nature of the crystals with dis- 
tance is predicted by chemical equilibrium and radial mix- 
ing models (Gail, 2004). The innermost disk is expected 
to contain mostly forsterite, while at larger distance from 
the star a conversion from forsterite to enstatite takes place. 
Thermal anneaUng produces both olivines and pyroxenes, 
with a ratio depending on the stochiometry of the amor- 
phous material accreted from the parent molecular cloud. 
van Boekel et al. (2004a) also find evidence for a decrease 
of the average grain size with distance from the star, sug- 
gesting that grain aggregation has proceeded further in the 
innermost, densest disk regions. 



7.3 TW Hya Systems 

A number of objects with distinct signs of evolution of 
their inner disk has been identified in the last decade. These 
objects with "transitional" disks are characterized by spec- 
tral energy distributions (SEDs) with a flux deficit in the 
near-infrared relative to the median SED of Classical T 
Tauri stars in Taurus (a proxy for the expected emission 
from optically thick disks, see D'Alessio et al., 1999), while 
at longer wavelengths fluxes are comparable or sometimes 
higher than the median SED of Taurus. Ground-based near 
and mid-infrared photometry combined with IRAS fluxes 
allowed the identification of objects with these properties 
in Taurus (Marsh andMahoney, 1992; Jensen andMathieu, 
1997; Bergin et al, 2004) and in the TW Hya association 
(Jayawardhana et al., 1999; Calve t et al., 2002). How- 
ever, instruments on board Spitzer are providing a much 
better view of the characteristics and frequencies of these 
objects (Uchida et al., 2004; Forrest et al., 2004; Muze- 
rolle et al, 2004; Calvet et al, 2005; Sicilia-Aguilar et al., 
2006). The inner disk clearing may be related to photo- 
evaporation of the outer disk by UV radiation (Clarke et al., 

2001) . However, transitional disks have now been found in 
brown dwarfs (MuzeroUe et al, 2006), for which the UV 
supply of energy is neghgible. An alternative explanation is 
that a planet has formed and carved out a gap in the disk; hy- 
drodynamical simulations seem to support this view (Rice 
et al, 2003; Quillen et al, 2004). 

The 10 Myr old TW Hya was the first of these objects 
analyzed in detail. The peculiar SED of this object (Cal- 
vet et al, 2002; Uchida et al, 2004), can be understood by 
an optically thick disk truncated at ~ 4 AU; the wall at the 
edge of this disk is illuminated directly by stellar radiation, 
producing the fast rise of emission at wavelengths > 10/xm. 
The region encircled by this wall is not empty. For one 
thing, it has gas, because the disk is still accreting mass 
onto the star (MuzeroUe et al, 2000). In addition, this in- 
ner region contains ~ 0.5 lunar masses of micron-size parti- 
cles, responsible for the small near-infrared flux excess over 
photospheric emission. This dust is also responsible for the 
strong siUcate feature at 10 /zm, which allows us to exam- 
ine the conditions of the dust in the inner disk (Calvet et al, 

2002) . The profile of the 10 /xm feature, integrated over the 
disk, is typical of amorphous siUcates with sizes of 2 /xm 
(Uchida et al, 2004), with a very low content of crystalline 
sihcates, < 2 % fraction by mass (Sargent et al, 2005). 
However, interferometric observations with MIDI resolve 
the feature from the inner disk (within ^1 AU; see Leinert 
et al. 2006 and the chapter by Millan-Gabet et al), show- 
ing clearly the 11.3 fjm peak typical of forsterite, while the 
uncorrelated profile is very similar to that measured with 
single-dish instruments. 

Sargent et al. (2006) have analyzed the dust content of a 
number of transitional disks, in addition to TW Hya. They 
find, averaged over the whole disk, a very low content of 
crystalline siUcates. In contrast, there is a large spread of 
the crystalline silicate fraction (from 0.3% to 20% in mass) 
in their sample of objects with optically thick inner disks. 
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These results support the idea that disks accumulate 
crystalline silicate grains in their innermost regions and that 
radial mixing transports outward some fraction of this mate- 
rial, with an efficiency that decreases with radius. In transi- 
tional disks, a large fraction of the region within ~few AU 
is cleared out, leaving only material with low crystalhne 
content. Spatially resolved observations of the innermost 
regions are needed to detect the small amount of remaining 
grains, with their higher crystalhnity. 

8. SUMMARY AND OPEN PROBLEMS 

8.1 Evidence of grain growth 

The results we have discussed come from a variety of 
techniques, including optical, near and mid-infrared imag- 
ing, mid-infrared spectrometry and milhmeter interferom- 
etry. All the data show clear evidence that grains in pro- 
toplanetary disks differ significantly from grains in the dif- 
fuse ISM and in molecular clouds. They are much larger 
on average, and there is increasing support for a scenario 
in which grains grow and sediment toward the disk mid- 
plane. On the surface of the disk, within few tens AU from 
the star, we know that in most cases silicates have grown 
to sizes of few microns, although much smaller particles 
(PAHs) may be present as well, and there are hints of ver- 
tical stratification (at several scale heigths from the mid- 
plane) of these micron-size grains from studies of scattered 
hght and silhouette disks. We have evidence from millime- 
ter interferometry that, if we consider not only grains on 
the disk surface but the bulk of the dust mass, in many disks 
grain growth has not stopped at micron-sizes, but there is 
a dominant population of "sand and pebbles", milhmeter 
and centimeter grains. In addition to vertical gradients of 
grain sizes, there is also evidence that radial gradients of 
grain properties exist, as shown by mid-IR interferometry, 
absorption mid-IR spectra of silhouette disks, and millime- 
ter images at different wavelengths, as discussed before. 
Grains in the very outer disk seem to be less processed than 
grains closer to the star. 

One important aspect to keep in mind is that not all the 
objects show evidence of processed grains: in some stars, 
for example, the 10 /im silicate feature has a shape very 
similar to that of the small silicates in the ISM; some disks 
have a rather steep dependence of the millimeter flux on 
wavelength, again typical of the small ISM grains. Some 
objects seem to have small sihcates on the disk surface, 
while very large grains are imphed by the milhmeter spec- 
tral shape. An impressive case is that of UX Ori, which 
has a very peaked and strong 10 /^m emission feature, typ- 
ical of the small ISM amorphous silicates (Fig. 2) and an 
extremely flat millimeter spectrum, which may indicate the 
presence of very large grains in the midplane. 

One should also remember that the very nature of pro- 
toplanetary disks (i.e., high optical depth, strong tempera- 
ture radial gradient, etc.) prevent a complete census of the 
dust population, so that deriving from the observations the 



global properties of the sohds in any given disk (such as 
the density of grains of different size and composition as 
function of radius and altitude) remains for the moment im- 
possible, and we have to do with a piecemeal picture. 

8.2 Dependence on stellar properties 

One important step forward since PPIV has been the ca- 
pabihty to study disks around stars of very different prop- 
erties. The advent of 10-meter class telescopes on the 
ground, the success of Spitzer and the improvement of the 
millimeter interferometers, including the VLA, have given 
us access to objects of increasingly lower mass. We have 
now millimeter spatially resolved data for TTS, and mid-IR 
spectra for very low mass stars and brown dwarfs, which al- 
low us to investigate if there is a dependence of grain evolu- 
tion on stellar and disk parameters. It may be worth noting 
that, within the millimeter-observed sample, which includes 
TTS and HAe stars, the stellar mass varies by a factor 5 
and the luminosity by a factor 250; assuming typical scal- 
ing laws (Muzerolle et al. 2003; Natta et al, 2004b, 2006), 
the accretion rate in the disk is likely to vary by at least a 
factor 25. The range of stellar properties are even larger 
for the mid-infrared sample, which includes now several 
brown dwarfs and covers the interval between ~ 0.04 and 
2-3 M0 in mass, ~ 0.01 and 50 L© in luminosity, ~ 10~^° 
and lO^*" M,-:/yr in accretion rate. Over this large range 
of physical conditions, we do not detect systematic differ- 
ences, i.e., grain properties do not seem to correlate with 
any star or disk parameter. At any given mass, there are 
large variations of grain properties, although objects with 
unprocessed grains seem to be rare. 

It seems that disks around all kind of stars have the po- 
tential of processing grains, to a degree that varies for rea- 
sons not understood so far. 

8.3 Time dependence 

The time dependence of grain growth is a crucial infor- 
mation for models, and one that we would like to derive 
from the observations with how much detail as possible. 
We should stress, first of all, that there is not, at present, 
anything close to an unbiased sample which can be used 
to study grain evolution. All the observing techniques dis- 
cussed in this review have been apphed to selected objects, 
known to be suitable for detection and analysis. In addi- 
tion, ages of individual pre-main sequence stars (our best 
clock so far) are measured from the location on the HR dia- 
gram and have large uncertainties, especially for HAe stars 
and brown dwarfs. Even more importantly, the objects for 
which we have observations vary in age by about a factor 
of ten only (from ^ 1 to ^ 10 Myr), and individual scatter 
may dominate over an underlying time dependence. 

The data available so far show no evidence of a correla- 
tion of the grain properties (at any wavelength) with time: 
we find very processed grains in the disks of the youngest 
as of the oldest pre-main sequence stars. One possible ex- 
planation is that grains are processed efficiently in the very 
early stages of disk evolution, when the system star-ndisk 
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is still embedded and accreting actively from the natal core 
(see, e.g., the theoretical models of dust dynamics and evo- 
lution during the formation of a protostellar disk by Suttner 
and Yorke, 2001). Then, the large star-to-star variation seen 
in, e.g., the HAe stars would be due to differences in pro- 
cessing during this initial, very active phase. However, we 
cannot exclude that further modifications also take place in 
the ~ 10 Myr in which the disk conlinues to exist. 

Some information can be obtained by studying grain 
properties in disks in Class I sources, which are in an ear- 
lier evolutionary stage than the disks around optically vis- 
ible objects discussed in this review. In Class I objects, 
the silicate features are seen in general in absorption, and 
show profiles typical of small, unprocessed amorphous sil- 
icates. It is likely, however, that the absorption is domi- 
nated by grains in the surrounding envelope, for which we 
expect very little growth. Only if the balance between the 
disk emission and the intervening absorption in the core is 
favorable, can we have a glimpse at the properties of the 
disk silicates, and in few cases there is evidence that the 
disk sihcates are much larger than those in the core (e.g., 
Kessler-Silacci et al. 2005). Also, the young protostellar 
binary SVS20 in Serpens shows evidence of the 11.3 /im 
forsterite feature in its spectrum (Ciardi et al, 2005). For 
the moment, there are only a few good quality spectra, but 
we can expect major progress soon, as Spitzer data will be- 
come available. One would also hke to know if Class I disks 
have milUmeter spectral energy distributions as flat as those 
of Class 11 objects. Unfortunately, using millimeter inter- 
ferometry to measure the spectral index of the emission of 
disks embedded in Class I cores has proved so far very diffi- 
cult, as the core emission is substantial even at the smallest 
physical scales one can probe. The extremely flat spectra of 
the few objects studied so far come from unresolved cen- 
tral condensations, and cannot be interpreted as evidence of 
very large grains in disks until the emission can be spatially 
resolved (Hogerheijde etal, 1998, 1999). 

8.4 Comparison with grain evolution models 

Models of grain growth by colUsional coagulation and 
sedimentation (e.g., DuUemond and Dominik, 2005; Wei- 
denschilling and Cuzzi, 1993, and references therein) pre- 
dict that grain growth will occur on very short time scales. 
This, at first sight, is in agreement with the idea that grain 
properties change already in the early evolutionary phases, 
when disks are actively accreting matter from the parental 
core. 

However, models also predict that the growth will not 
stop at micron or even at centimeter sizes, but will proceed 
very quickly (compared to the pre-main sequence stellar life 
times) to form planetesimals. This is not consistent with 
the observations, which show us disks around stars as old 
as 10^ years with grains that, albeit very large, are far from 
being planetesimals. When fragmentation of the conglom- 
erates is included, in many cases the outcome is a bimodal 
distribution for the solids, with most of the mass in planetes- 
imals and a tail of smaller fragments, which can be as large 



as centimeters and meters, but contains only a tiny fraction 
of the mass in planetesimals. If these smaller bodies pro- 
vide sufficient optical depth, the disk properties may look 
similar to what is observed. However, we think this is un- 
likely to be the case. The millimeter observations described 
in Section 6 tell us not only that grains need to have very big 
sizes but also the mass of these grains, which turns out to be 
very large. Using, e.g., the same grain models as in Fig. 5, 
we derive dust masses, for the disks with evidence of very 
large grains, between 10~^ and 10~^ M©; assuming the 
standard gas-to-dust ratio of 100 (which is probably close 
to true when disks formed), this implies disk masses be- 
tween ~ 0.1 and 1 M0, i.e., ratios of the disk to star masses 
approaching unity. If a large fraction of the original dust is 
not in the observed population of millimeter and centimeter 
grains, but in planetesimals, the disk mass should be even 
higher, well above the gravitational instability limit. Al- 
though we cannot exclude that a small fraction of the orig- 
inal solid mass is in planetesimals, we think it is unlikely 
that they have collected most of it. 

The survival of a large mass of milhmeter and 
centimeter-size grains over a time scale of several million 
years, as observed in many disks, is challenging the models 
on several grounds. It suggests that in many cases the pro- 
cess of planetesimal formation is much less efficient than 
predicted by theoretical models. If disks are forming plan- 
etesimals as fast as predicted, then this can only happen at 
the end of a long quiescent phase, where growth is limited 
to sizes of few centimeter at most; otherwise, the process 
has to be very slow, involving a modest fraction of the dust 
mass for most of the pre-main sequence life of the star. The 
"inefficiency" may result from fragmentation of the larger 
bodies. This in-situ production of the grains we observe 
may take care of the other severe difficulty one has, namely 
that millimeter and centimeter solid bodies (in a gas-rich 
disk) migrate toward the star on very short timescales. This 
is a very open field, where progress can be made only by 
combining together theoretical studies, laboratory experi- 
ments and observations. 

It is possible that the disks we can study with current 
techniques are those that will never form planets. Disks 
where planets form may evolve indeed very quickly, and 
could be below present day detection limits. This, how- 
ever, seems unlikely, as we have objects like TW Hya and 
HD 100546, which have evidence of very large grains in 
their outer disk, and an inner gap possibly due to the action 
of a large planet. Similar "transitional" objects are found 
by Spitzer in increasing number, and a detailed character- 
ization of their millimeter properties will be important. In 
any case, the properties and evolution of disks that we have 
studied so far and their dust content provide the only obser- 
vational constraints to theoretical models of grain evolution 
and planet formation. 
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8.5 Crystalline silicates and mineralogy 

The ISO discovery of large fractional amounts of crys- 
talline silicates in disks around stars of all masses has been 
a surprise. They are enhanced in the innermost regions 
of disks, as shown by interferometric observations of HAe 
stars and by the spectra of transitional disks, but they are 
also present, in variable degree, further out. Fractional 
masses of crystalline sihcates, integrated over the whole 
disk emitting regions, have been derived for many objects, 
HAe, TTS and BDs, with values ranging from zero to about 
60%. These estimates refer to grains on the disk surface, 
typically within 10-20 AU from the star for the HAe stars, 
and 0.1-0.2 AU for BDs. One should note that these values 
are extremely model dependent, and refer to a limited range 
of grain sizes only; they should be taken with the greatest 
care. We can expect great progress in the immediate fu- 
ture, as spectra in the 15-45 yum spectral region, which can 
detect colder (and larger) crystals, are becoming available. 
Further insight can also come from improved characteriza- 
tion of the mineralogy in disks, e.g., measuring the ratio of 
species such as crystalline pyroxene, forsterite, sihca etc. 
Self-consistent radiation transfer models in disks which can 
include a large number of dust species exist, and can be used 
to provide more rehable values for the relative abundances 
of the various species observed on disk surfaces. 

The formation of crystalline silicates requires energetic 
processes, which have to be efficient not only in HAe stars, 
where crystalline silicates were first discovered, but also in 
the much less luminous TTS and brown dwarfs. If crystal- 
hzation is restricted to the very inner disks, as in the Gail 
(2004) models, then strong radial mixing must occur, trans- 
porting material outwordly. Radial and meridian drifts can 
have an important role in the growth of grains, which need 
to be addressed in the future. 

8.6 Final remarks 

The observations we have discussed in this chapter show 
clearly that very few disks (if any), contain unprocessed 
grains, i.e., with properties similar to those of grains in the 
ISM. In general, dust has been largely processed in all ob- 
jects. 

However, the degree and the result of these changes vary 
hugely from object to object, as we find large differences 
in the grain size, composition and allotropic properties. At 
present, we have not been able to identify any correlation 
between the grains and other properties of the central star 
(such as its mass and luminosity) and of the disk (for ex- 
ample, mass, accretion rate, etc.), nor with the age of the 
system. This last point is particularly distressing, since 
planet formation theories require to understand if, when and 
how grains change with time, or, in other words, how they 
"evolve". 

The study of pre-main sequence stars has shown us that, 
no matter which aspect of their rich phenomenology we are 
interested in, there is a large scatter between individuals, 
which can hide underlying trends. In discussing grain prop- 
erties, it is possible that the objects observed so far do not 



sample the right range of parameters, age in particular, or 
that the available samples are still too small, or restricted to 
a too narrow range in the parameter space. We see the noise, 
and cannot identify trends which remain, if present, hidden. 
If this is really the case, we can expect major advances in 
the near future, as new space and ground-based facihties 
will make possible to study much larger samples of disks, 
in different star-forming regions, surrounding objects dis- 
tributed over a broader range of mass, age and multiphcity. 
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